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The criteria for validity of adiabaticity for nonlinear wave equations are considered within the 
context of atomic matter-waves tunneling from macroscopically populated optical standing-wave 
traps loaded from a Bose-Einstein condensate. We show that even when the optical standing wave 
is slowly turned on and the condensate behaves adiabatically during this turn-on, once the tunneling- 
time between wells in the optical lattice becomes longer than the nonlinear time-scale, adiabaticity 
breaks down and a significant spatially varying phase develops across the condensate wave function 
from well to well. This phase drastically affects the contrast of the fringe pattern in Josephson-effect 
interference experiments, and the condensate coherence properties in general. 



Introduction. — Recent matter-wave interference studies with Bose-Einstein condensed atoms have shown that 
one can observe phase-dependent dynamics in dilute neutral-atom systems, in complete analogy with those observed 
in Josephson-effect experiments |J. In some of the experiments reported in Refs. |3|,[|, the external potential (applied 
by optical standing waves to the initial ground-state Bose-Einstein condensate (BEC)) is slowly turned on. One 
might assume, based upon the Adiabatic Theorem of Quantum Mechanics , that the system should remain in an 
eigenstate, even though the nature of the eigenstate evolves in time. Moreover, adiabaticity criteria for nonlinear wave 
equations have been studied in connection with soliton dynamics in nearly integrable systems, and one expects on 
the basis of these studies that adiabaticity should be maintained with slow enough variation of the parameters of the 
system Q . A spatially varying phase of the condensate wave function is a manifestation of non-adiabatic dynamics 
of a BEC, in the sense that the adiabatic eigenstate (the ground nonlinear eigenstate of the Gross-Pitaevskii equation 
(GPE) calculated using the potential at any instant of time) can be taken to be real (it does not have a spatially 
depenendent phase). Hence one concludes that, a spatially varying phase should not develop across the condensate 
in Refs. (^^] when the optical potential is turned on very slowly. Here we present calculations corresponding to 
conditions similar to those reported in showing that the BEC remains adiabatic as the light potential is turned 
on very slowly, and the phase is constant across the condensate. The dynamics are indeed adiabatic despite the fact 
that nonlinear dynamics precludes the possibility of a superposition principle, which is used heavily in deriving the 
Adiabatic Theorem. We further show that, as the light potential is slowly increased in strength to the point where 
the wave packets in the individual optically-induced wells become strongly separated, and the time-scale for tunneling 
becomes long compared to the nonlinear time-scale, adiabaticity is destroyed and a large spatially varying phase 
develops across the wave function. This large inhomogencous phase, obtained within a mean-field approximation, is 
due to the tunneling dynamics (induced by the kinetic energy of the BEC in the lattice potential) not being able to 
implement the equilibration of the phase across the BEC on the nonlinear time-scale. 

The process of dynamic splitting of a condensate by an external time- varying potential was studied using two mode 
models [|l(| which show that the mean-field approximation based upon the GPE will not give a good description of the 
splitting process due to the slow rise of a potential barrier which cuts off tunneling processes. Our results show that, 
even within a mean-field approximation, a sudden inhomogeneous phase build-up begins to develop when the tunneling 
time becomes comparable to or larger than the nonlinear time-scale, and the process of optical lattice turn-on ceases 
to be adiabatic when this happens. A transition occurs when these two time-scales become comparable, already 
within a mean-field GPE approach. The spatially varying phase adversely affects the fringe contrast in interference 
experiments performed on the BEC wavepackets. 

Adiabaticity in nonlinear systems can be studied in various regimes. Denoting tad as the quantum-mechanical 
linear adiabatic time-scale determined in terms of the difference of the (linear) energy eigenvalues and tjvl the 
nonlinear time-scale [ pd| (see below), the simplest regime is one in which the the duration of the dynamical process 
being studied, T, satisfies the condition, tad <C T <C tml- In this case, adiabaticity is insured by the Adiabatic 
Theorem due to the first inequality, and nonlinearity cannot play a significant role in the dynamics due to the second 
inequality. Hence, the dynamics must be adiabatic. The regime in which the experiments of Refs. are carried out 
satisfy the condition tad, Tnx <T, and the nonlinearity does play an important role in the dynamics. Nevertheless, 
as described below, the dynamics are indeed adiabatic, until the strength of the optical lattice is so large that the 
wave packets in the individual wells become almost completely separated. 
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We consider Bose- Einstein condensed 87 Rb atoms in the \F — 2,Mp — 2) hyperfine state confined in an array of 
optical traps in a gravitational field. The atoms are trapped at the antinodes of a vertically oriented red-detuned 
optical standing wave, which are separated by Az = A/2, where A = 840 nm is the wavelength of light used to confine 
the atoms. The depth of the optical potential is proportional to the intensity of the light; the intensity is initially zero 
and linearly increases with time. The initial BEC is cigar shaped with 10 4 atoms in a static harmonic trap potential 
V ho (r) = ^z 2 + ^^-{x 2 + y 2 ) with frequencies tu z = 2tt x 19 Hz, and lu = 2n x 33 Hz (uj = {uj^WywJ) 1 ^ with 
u> x = uj y = uj xy ) B. The light-potential experienced by the atoms in the BEC is given by Vl(z, t) — Vb(t)(l+cos(2fc£z)) 

where the well-depth, Vo(t) = 125 t [En], varies linearly with time. The recoil energy Er = 2n ^ is the kinetic energy 
gained by an atom absorbing a photon from the optical lattice, where m is the atomic mass, and the photon wavevector 
is fci = 2tt/X. The rate of increase of the light-potential, (125 [Er/ s]), is sufficiently slow, as shown below, that the 
dynamics of the BEC is adiabatic for much of the turn-on. After some time, the harmonic potential and the light- 
potential are switched-off (dropped) releasing the atoms to fall under the influence of gravity. In the experiments 
reported in Ref. 0], the atoms are held in the optical lattice for a short time (2.5 ms) after switching off the harmonic 
potential, allowing the gravitational potential difference between wells to affect the phase difference between wells, 
and absorption images are taken 8 ms after the optical potential is turned off and the atoms begin to free-fall. 

Theoretical Formulation — The mean-field dynamics can be determined in terms of the time-dependent CPE, 
ih^j^ 1 = (p 2 /2m + V(r,t) + N U \if)] 2 )tp, where V(r,t) = V ho (r) - mgz + V L (z,t), U = is the atom-atom 

interaction strength that is proportional to the s-wave scattering length oq , and Nq is the total number of atoms. We 
solve the time-dependent GPE using a split-operator fast Fourier transform method to propagate an initial state of 
the BEC in time in the presence of the harmonic potential, gravity and the time dependent optical lattice; the initial 
state is determined by propagating in imaginary time with vanishing optical lattice potential p"j] ]. Due to the large 
number of grid points necessary in the lattice direction (z), and the large number of time steps necessary to propagate 
the GPE to completion of the dynamics, we found it necessary to convert the 3D GPE into an effective ID GPE with 
similar dynamics. This is carried out using the following procedure. The ID GPE is written in terms of characteristic 
time-scales tr>F for diffraction and tj^L for the nonlinear interaction, in the following manner 
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Here, rxF = \/2/V (mui 2 ) is the Thomas-Fermi radius, /1 is the initial chemical potential, fi — \ ( 15 ^ jV ) (mu) 2 ) 3 / 5 , 
tjjp is the diffraction time tup — 2mrj, F /h and [ip m \ is the maximum initial magnitude of the wave function defined 
in terms of the initial nonlinear time by ijvL = (/V^) 1 = (NUo\4> m \ 2 /h) 1 PI" ^ e all lengths in the GPE (|l|) 
in units of the Thomas- Fermi radius along the z axis, and we call this length unit r z = Ttf,z- The key to obtaining 
physically relevant dynamics using the ID GPE, is to use (1) uj\d = (^f - ) 1 ^ 2 ; an d (2) N±d = Nq(u)/ujid) 3 - These 
two equations insure that the Thomas-Fermi radius and the nonlinear time (or the chemical potential) remain as in 
3D world. We checked to confirm that this procedure gives the same wave function as the ID projection of the 3D 
GPE solution for our studies using shorter propagation times. 

Numerical Results — Starting from the BEC without any optical lattice present, we begin to increase the optical 
potential with the linear ramp Vo(t) mentioned earlier. Fig. |l| shows the magnitude and phase of the Gross-Pitaevskii 
wave function, ip(z, t) = \ip(z, t)\ exp(i#(z, t)) as a function of position in the optical lattice, z/r z , when the well-depth 
of the optical lattice is 10 and 15 Er. We have taken z — to be at the minimum of the harmonic potential plus the 
gravitational potential, i.e., we have shifted the definition of the center of the trap to the true center of the combined 
harmonic plus gravitational potential. About 17 wells of the optical potential are populated during the course of 
the dynamics. The magnitude \4>(z,t)\ remains normalized throughout the propagation; if we were to average out 
the oscillations in |^>(z, t)| we would obtain roughly the initial BEC \ip(z, 0)|. The phase 9(z,t) is nearly spatially 
independent, and the magnitudes \ip(z, t)\ in Fig. |l| are very nearly equal to the eigenstates of the GPE with the 
optical potential present at these times. The slight variation of the phases with position z indicate the small degree 
of the nonadiabaticity that results during the dynamics. The phases 9{z,t) vary as a function of time by a spatially 
independent constant due to the dynamics but this (z-independent) constant phase is not physically significant in 
the experiments of Refs. |M. The phase 6{z,t) at the time when the depth of the optical lattice is 10 Er has been 
shifted up in Fig. [l] by unity so that it could be plotted conveniently. 

As we increase the well-depth of the optical lattice further, a regime is reached in which the wave packets localized 
in the various wells become almost completely separated. When this occurs, the spatially varying phase of the wave 
function begins to grow significantly, and the spatially dependent variations increase as the well-depth of the optical 
lattice increases. Fig. || shows the magnitude \ip(z, t)\ and phase 0(z,t) of the GP wave function as a function of 
position in the optical lattice, z/r z , when the well-depth of the optical lattice is 25 and 40 Er. The magnitude of the 
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wave function is almost totally within the region z/r z E [—1-3, 1.3]. The wave packets in the optical wells are almost 
fully separated. The spatially dependent variations in the phase are less pronounced in the center of the trap where 
the variation in the peak densities of the wave packets from well to well are smaller; by 40 Er the variation in phase 
across the condensate is large compared with 1 radian. The jump in the phase for 40 Er at zjr z « ±1.1 is artificial 
and due to the continuation of the inverse trigonometric function used to calculate 9(z,t) from ip(z,t). The spatial 
phase variation over the condensate for 50 Er is much greater than that for 40 Er; it varies over more than 5w. 

Fig. ^ shows the magnitude of the Fourier transform of the Gross-Pitaevskii wave function ip(k) versus k when the 
well-depth of the optical lattice is 25, 40 and 50 Er. The wave function has amplitude around k = 0, ±2&£, although 
components around k — ±4kL are also clearly visible. The width of the wavepackets increase with increasing well- 
depth, particularly as the spatially varying phase across the condensate becomes significant . By 50 Er, the width of 
Fourier components of the wave packets are a good fraction of 2fc^. 

Adiabaticity is maintained throughout the course of the dynamics until the wave packets become well separated and 
the tunneling time becomes comparable to or larger than the nonlinear time-scale. Adiabaticity of soliton solutions 
of nonlinear wave equations have been extensively studied for slowly varying external conditions (for a review, see 
Ref. H). It has been shown that if the time-scale T of the variation of an external parameter is slow compared the 
instantaneous nonlinear eigenvalue LUo(t) of the nonlinear (time-independent) equation at time t, u>o(t)T / '(2tt) 3> 1, 
the dynamics can be adiabatic. Here we see that an additional condition is required; adiabaticity breaks down when 
the wave packets become well separated and the tunneling time becomes longer than the nonlinear time. 

It is easy to make a rough estimate of the tunneling probability from well to well using the semiclassical approxi- 
mation ||, 

/ W(2fci) \ / [2^7 r /[2kL) , \ / , \ 

P(t) = cxpi-J k(z,t)dz\ ^expl-y^- J (V L (z,t))^ 2 dz\ = cxp (-2y/4mV (t)/(hk L )) . (2) 

The tunneling rate is then given by R(t) = u) v {t) x P(t), where the vibrational frequency u> v (t) can be approx- 
imated by uj v (t) w \J Ak'jVoif) I 'm, since the expansion of the optical potential about a minimum in the poten- 
tial yields Vh{z,t) w (m/2)(uj v (t)) 2 z 2 = (m/2)(Ak 2 j Vo(t)/m)z 2 . Hence, the time-dependent tunneling rate is 

R(t) — y/8(Vo(t)/ER) exp (^—^8(Vo(t)/ER)^j [Er/H]. The tunneling rate R(t) can be compared with the time- 
dependent inverse nonlinear time-scale, (tNLit)) -1 = NUo\ifj m (t)\ 2 /h = fi(t)/h and the time-dependent diffraction 
time-scale, (i_DF(i)) -1 — 4j^fr ■ Once the tunneling time, R~ l , becomes long compared to ijvL (and tjjp, which is 
typically the smallest of these time-scales), a spatially independent phase can not be maintained across the conden- 
sate by the action of the kinetic energy operator, and a spatially dependent phase develops. In our calculations, this 
happens beyond a lattice well-depth of about 35 Er. 

If we take our calculated wave packets (shown in Figs. [I] and |2|) and propagate for an additional 2.5 ms upon 
switching off the magnetic field, but leaving the optical lattice potential in place, as described in Ref. [0, our wave 
packets hardly change. However, upon switching off the optical potential and propagating the wave packets for 8 
ms, the wave packets spread by diffusion very significantly as they fall in gravity. The resulting wave packets are 
considerably different upon using the 25 and 50 Er results. The additional Fourier components of the wave packet 
in the 50 Er case significantly wash out interference patterns in the density profile for this case, consistent with the 
measurement in Note that we are not asserting that squeezing is absent in the experiments of Ref. 0]; rather, 
that an improved model of the above-the-mean field effects that accounts for the spatially varying phase is necessary 
to quantitatively compare with the experiments. 

Conclusions — We have seen that the spatially varying phase of the solution to the Gross-Pitaevskii equation 
becomes large when the tunneling time between wells in the optical lattice becomes comparable to or larger than 
the nonlinear time-scale. Adiabaticity then breaks down and the instantaneous nonlinear eigenvectors to the time- 
independent Gross-Pitaevskii equation do not have the character of the dynamical solution. This nonadiabaticity takes 
place even when the variation of the optical potential is sufficiently slow that adiabaticity is otherwise assured. It 
remains to determine the above-mean-field corrections to this picture using an approach that incorporates an accurate 
form for the spatially varying complex order parameter. 

Finally, we note that these results have implications regarding the loading of an optical lattice with atoms for 
quantum computing, using a Bose- Einstein condensate (BEC) source and laser beams that are slowly turned on [ jT3"| . 
From our studies we conclude that the laser fields must be controlled so that the density of the atoms is sufficiently 
low (the nonlinear term must be negligible) before tunneling between the wells is cut off by the full depth of the optical 
potential. Otherwise, phase variations from well to well will deleteriously affect the resulting optical lattice state. 



3 



Useful conversations with Boris Malomed and Paul Julienne are gratefully acknowledged. This work was supported 
in part by grants from the U.S. -Israel Binational Science Foundation (grant No. 98-421), Jerusalem, Israel, the 
Israel Science Foundation (grant No. 212/01), the Israel MOD Research and Technology Unit and the Polish KBN 
2/P03/B07819 



[1] M. R. Andrews et al, Science 275, 637 (1997). 

[2] E. W. Hagley, et al, Phys. Rev. Lett. 83, 3112 (1999); M. Trippenbach, et al, J. Phys. B 33, 47-54 (2000). 

[3] B. P. Anderson and M. A. Kasevich, Science 282, 1686 (1998). 

[4] C. Orzel et al, Science 291, 2386 (2001). 

[5] F. S. Cataliotti et al, Science 293, 843 (2001). 

[6] A. Barone and G. Paterno, Physics and Applications of the Josephson Effect (Wiley, New York, 1982). 
[7] A. Messiah, Quantum Mechanics, Vol. II, Chp. 17, (N. Holland, Amsterdam, 1975). 
[8] C. Cohen- Tannoudji, B. Diu and F. Laloe, Quantum Mechanics, (John Wiley, NY, 1977). 
[9] Y.S. Kivshar and B.A. Malomed, Rev. Mod. Phys. 61, 763 (1989). 
[10] C. Menotti, J.R. Anglin, J.I. Cirac and P. Zoller, Phys. Rev. A63, 023601, 2001; A. J. Leggett and F. Sols, Phys. Rev. 
Lett. 81, 1344 (1998); R.W. Spekkens and J.E. Sipe, Phys. Rev. A 59, 3868 (1999); J. Javanainen and MYu. Ivanov, 
Phys. Rev. A 60, 2351 (1999); A. Vardi and J.R. Anglin, Phys. Rev. Lett. 86, 568 (2001). 
[11] M. Trippenbach, Y. B. Band, and P. S. Julienne, Phys. Rev. A62, 023608 (2000). 
[12] M. Trippenbach, Y. B. Band, and P. S. Julienne, Optics Express 3, 530 (1998). 

[13] O.K. Brennen, CM. Caves, P.S. Jessen, and I.H. Deutsch, Phys. Rev. Lett. 82, 1060 (1999); D. Jaksch, et al, Phys. Rev. 
Lett. 82, 1975 (1999). 




FIG. 1. Magnitude and phase of the Gross-Pitaevskii wave function as a function of position in the optical lattice. The wave 
function is shown when the well-depth of the optical lattice are 10 and 15 Er. 
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FIG. 2. Magnitude and phase of the Gross-Pitaevskii wave function as a function of position in the optical lattice. The wave 
function is shown when the well-depth of the optical lattice are 25 and 40 Er. 
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FIG. 3. Gross-Pitaevskii wave function in momentum space, ip(k) versus k when the well-depth of the optical lattice are 25, 
40 and 50 E R . 
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